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We describe the application of differential reduction algorithms for Feynman Diagram 
calculation. We illustrate the procedure in the context of the generalized hypergeomet- 
ric functions P +iF p , and give an example for a type of g-loop bubble diagram. 

1 Introduction 

A Feynman diagram can be understood mathematically as a linear combination of Horn- type 
hypcrgcometric functions of several variables (see Ref . [5] ) : 

3 

where a Horn-type hypergeometric function has the structure 

*/ V fUf=l^(E r a=1^3ama+-f 3 )\ 

$(7; a; x) = > — j \x" 1 ■ ■ ■ x™ r , (2) 

mi,m a ,... ,mr=0 \ Uk=l T (Eb=l v k b m b + a k ) ) 

with the arguments Xj being, in general, rational functions (typically, simple ratios) of 
kinematic invariants of the original Feynman diagram, and the parameters {7^} and {ofcl 
being linear combinations of the exponents of propagators and the dimension of spacc-time|f| 
The 7j and a k are called upper and lower parameters, respectively. 

These statements follow from the multiple Mcllin-Barncs representation for a dimension- 
ally regularized Feynman diagram (see Ref. [?]), and the assumption that there is a region 
of variables where every term in the linear combination ([1]) is convergent. The Horn- type 
structure permits the hypergeometric functions appearing in |T]) to be reduced to a set of 
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a The presence of a nontrivial numerator in the Feynman diagram docs not affect this conclusion [3]. 
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basis functions with parameters differing from the original ones by integer shifts: 



P (x)$(j + ka + s;x)= £ P ri ,.. ,r p (x)D\ n - ■ ■ ■ D^H(f, a; x) , (3) 

mi,-- - ,m r —0 

where = fgf - ) denotes a partial derivative and P rit ... tTp (x) are rational functions 

[5]. These shifts may be implemented by constructing a set of four differential operators 
U^ c , L^ a which respectively change y c , cr c by ±1: 7 C — *■ 7 C ± 1 or er c — > ct c ± 1. These 
basic operators are called the step-up and step-down operators for the upper and lower 
parameters. A procedure of applying step-up and step-down operators to reduce the original 
hypergcomctric function to a basis set is called a differential reduction. In the case when 
some of the variables are equal to one another, Xi = Xj,i ^ j, or belong to the surface of 
singularities Q = {x\Po(x) = 0}, it is necessary to define a limiting procedure for ©. 

We will illustrate our approach by considering the reduction of a particular generalized 
hypcrgeometric function in section [2j and by applying it to a particular class of Feynman 
diagrams in section [3] Section [4] describes the reduction at a singular surface using the e 
expansion and hyper logarithms. 



2 Generalized hypergeometric function of one variable 

In this section, we will show how differential reduction may be applied to a generalized 
hypergeometric function of one variable. Let us recall that the generalized hypergeometric 
function p i 7 p _i(a; b; z) may be defined in a neighborhood of z = by the series 



F{a- b-z) 



pFp-i 



fe=0 



(4) 



where (a)k = T(a + k)/T(a) is called a Pochhammer symbol. The lists a = {a\, ■ ■ ■ , a p ) 
and b = (pi, • ■ ■ , b q ) are called upper and lower parameters of the hypergeometric function, 
respectively. The hypergeometric function satisfies a differential equation 



~j{a,b) (^ p F p -i(a;b; z) 



IK 



p-i 

i=l 



p Fp-i(a;b;z) = 



(5) 



where L(a,b) is a differential operator and 9 = z-^ . 

Constructing a reduction scheme requires a set of step-up and step-down differential 
operators for both the upper and lower parameters. In this case, the universal step-up 
(step-down) operators for the upper (lower) parameters have a very simple form0 



J b 3 



6,-1 



and the inverse operators [/- , and Lt can be constructed in accordance with Takayama's 



algorithm, j^j (See also Ref. [2].) 



b See Ref. [5], or Eq. (2.1), (2.2) in Ref. [2], or Ref. [g] for details. 
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The differential reduction algorithm takes the form of a product of several differential 

±. 

b ' 



step- up/step-down operators U^,L^ 



D+1 F p (a + m;b + n;z) = (C/j)^™' ( L f) ' ' P +i F P (a; b; z) , 



(6) 



so that the maximal power of 9 in this expression is r = m, + ■ rij . Since the hypcrgc- 
ometric function p -\.\F p (a; b; z) satisfies a differential equation (J3J) of order p, it is possible to 
express all terms containing powers of 9 k with k > p in terms of 9^ p +\F p {a; b; z) with j < p, 
multiplied by coefficients that are rational functions of the parameters and the argument 
z. In this way, any function p +\F p (a + m;b + k; z) may be expressed in terms of a basic 
function p+ iF p (a; b; z) and its first p derivatives: 



S(a,i, bj, z) p+1 F p (a + m; b + k; z) = 



(7) 



Ri(di, b 3 ,z)9 p + R 2 {ai, bj, z)9 v 1 H h R P {a,i, bj,z)6 + Rp + i(ai,bj,z) > p+ iF p (a; b; z) , 



where to, k are lists of integers and S and T t are polynomials in the parameters {a^}, {bj} 
and z. For some special sets of parameters, the result of the reduction ([7]) takes a simple 
form. In particular, when one of the upper parameters is an integer, 

S{a i7 bj, z) p+ iFp(j, a+m; b+k; z) = 

Ri(ai, bj,z) + i^R 2 (ai, bj, z)9 p ^ 1 H h Rp{ai, bj,z)9 + R p+ i(ai, bj,z]^ p+ iF p (a; b; z) . 

For further details see Ref. [7]. 




3 Example of reduction 

As an example of how the differential reduction applied) 
to a particular type of Feynman diagram, let us consider 
the g-loop bubble diagram Bf 122 QQ in Fig. [TJ with four 
massive lines (with masses to, M as indicated) and two Figure 1: q-loop bubble -B 112 2oo 
sets of massless subloops with r and x lines, respectively. 
It is defined as 

B 11220o( m2 ' M2 ' a l' a 2^1 I / 3 2,CTi, • • • ,a x ,p 1 ,--- ,p r ) 

d n (k 1 ---k q ) 



KY 1 ■ ■ ■ [k 2 r _ 1 ]P^[(k 1 + - ■ ■ + k r ^ 1 + P) 2 ]Pr 



+x 



}^[{K+---+k r+x +py -to 



21 a 2 



(8) 



[k 2 ]^ ■ ■ ■ [k 2 r+x ^ [k 2 - Mf! \kl - Mf^ 



c See also Ref. [8] for more examples. 
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where P = + k q and x + r = q — 2 , g > 3 The Mellin-Barnes representation for this 
diagram it0 



- B 112200( to2 '- M2 ' Q; 1: Q; 2,/3i,/?2,O'1, " ' ,(T X ,PW 



,Pr) 



r(9-i)-^-/3i-/3 2 



r(ai)r(a 2 )ro9 x )ro%)r (f 



nr(5-ft)lr(p-5(r-i)) 



.i =1 



r Sr-p 



y ds r (a-f s-s) r (A+^+p+a-f (g-l)-s) 

r(ai+s)r(a 2 +s)r(a 1 +a 2 -|+s)r(f+s)r(|(g-2)- i o-cr+s) 

r(a 1 + a 2 + 2s)r(§(a; + l)-CT + .s) 
r(/3 1 +p+g-a( g -2)- s )r( / 3 2 + p+ < 7-§( g -2)-a) 
r(/3i+/3 2 + 2 ( 7+2 / 9-n(q-2)-2s) 



(9) 



where we have introduced the notations p = J2a=i Pa , & = ELi a a • 

Closing the contour of integration on the left (see Ref. [9] for details), the result may 
be expressed as a sum of four hypergeometric functions of type rF e . The general result 
is too long to reproduce here. In the case that the exponents of the propagators are all 
integers, the hypergeometric functions appearing in the result are reducible to the following 
hypergeometric functions and derivatives thereofEl 



{1,0,0 2 ,0 3 } x 4 F 3 
{l,0,0 2 ,0 3 }x 4 F 3 
{l,0,0 2 ,0 3 }x 5 F 4 



Ji-§(z-l),/ 2 -fz,/ 3 -i(z+l),| + /4+f(<?-z-2) 

i + / i-f' / 2-f(<Z-2),/ 3 -f(g-l),/4-f9 
h-^(q-x-l),I 6 -^(q-x-2),i+I 7 -^(q-l) 

i,i+/ 1 ,/3-f(g-i),/2-f(?-2),/ 4 -|(g-3) 

I 5 + %±+I 6 -%(q-2),I 7 -%(q-x-2)J 8 -%(q-x-3) 



(10) 



where Ik are arbitrary integers. In the last expression, some polynomials are also generated. 

Further simplification is possible when q, x take particular values. For example, for x = 
and x = 1, the first hypergeometric function is reducible to 



or to 



{1,0} X 3*2 



{l,0,0 2 }x 4 F 3 



i,/ 2 -f,±+/ 3 +f(g-2) 

I 4 + §(«-!) A+h 



z I (for x = 0) , 



%,h-n, 2 



-h + %{q-'i) 



z I (for x = 1) 



(11) 



(12) 



d We note that a smooth limit exists for <j — » and p — > 0. 

°Two of the original hyergeometric functions have a similar parameter structure and the same basis 
functions. 
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4 Reduction at a singular surface via the e expansion and hyper- 
logarithms 

In physical applications, the case of of equal masses m 2 = M 2 (a "single-scale" diagram) is 
of special interest^ For the diagram in Fig. |T]), this case corresponds to z = 1, which is a 
singular point for the differential reduction algorithm. The question is then how to find a 
smooth limit at this point. Let us recall [TT] that the hypergeometric function (H|) converges 
at z = 1 if Re (S6j — £dj) > 0. In this way, if the hypergeometric function on the l.h.s. of 
Eq. ([7]) is well-defined at z = 1, a smooth limit of the differential reduction exists. It is now 
a technical problem to rewrite the r.h.s. of Eq. ([7]) in terms of the variable x = 1 — z. It is 
well-known, however, that for p > 3, the hypergeometric function p F p _i is not expressible 
in terms of hypergeometric functions of the same type in the neighborhood of z = 1 (see 
Ref. [13]). 

One approach to this problem is to construct the all-order e-expansior@ in terms of 
functions which are defined for the entire range < z < 1. The problem is then solved at 
each order in e. The hyperlogarithms [TB] belong to that class. Unfortunately, at present, 
the necessary theorems on the all-order e-expansion are proven only for special sets of 
parameters [UJ [HI [19] . Let us recall some results from Ref. [18] , There are bases (sets 
of parameters) for which the all-order ^-expansion of a hypergeometric function has the 
form (see Ref. [18]) 

p F p - 1 (A + £a;B + eb;z) = C(A,a,B,b,z)J2 £ ' J2 c ? (a, 6)Li ? ( A*/, A^) , (13) 



J.k.s=l 



where X)i s * = 3 t ^ < ja, k < q, and q is integer number, the coefficients cg(a,b) are poly- 
nomials in the parameters {cij} and {b k }, X q is primitive g th -root of unity, Lij* ^X q A^£ 

are hyperlogarithms, \ q is short-hand for A^ 1 ~- 72 , Aj 2 ~ J3 , ■ ■ • , Xq 1 "^ 1 Jfc , £ is a variable re- 
lated algebraicallj0 to z, and C{A^a,B 1 b,z) is a polynomial. Thus, at each order in the 
e-cxpansion, only hyperlogarithms of a single weight arc generated. The definition of the 
hypcrlogarithm as an iterated integral over any rational function 



du. r h dt k _ 1 r 2 dt 



I(z;a k ,a k _ 1 , . . . ,a x ) 



where the dj are arbitrary numbers, together with the structure of the expansion (|13[) . 
show that the transformation z — > 1 — z results in functions of the same structure. In this 
way, we can construct the necessary limiting procedure without detailed knowledge about 



tk — a k Jo *fc-l ~ a fc-l JO h ~ fl l 

dt I(t;a k -i, •••;Oi) , 
o t — a k 



f See Ref. 1101 for the three-loop case (q = 3, x = 0, r = 0) and Ref. 1121 for the four-loop case (q = 4, x = 
l,r = 0). 

s See Ref. 1141 for another technique of evaluating hypergeometric functions at z = 1. 

i i / si 

For completeness, we note that £ can take the explicit forms £1,2,3 = z 1 , (1 — z) 1 , ( ) q where q is 
an integer. Under z — > 1 — z, these variables transform as £1,2 — * £2,1 and £3 — » ■ 
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the relationship between the functions p F p -\(a; b; z) and p F p -i(a;b;l — z). For practical 
applications to diagrams presently of interest, a few coefficients of the e-expansion suffice, 
and these have been implemented in several existing packages |20j . 
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